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ABSTRACT 

By combining test-particle and self-consistent techniques, we have developed a method to rapidly explore 
the parameter space of galactic encounters. Our method, implemented in an interactive graphics program', can 
be used to find the parameters required to reproduce the observed morphology and kinematics of interacting 
disk galaxies. We test this system on an artificial data-set of 36 equal-mass merging encounters, and show that 
it is usually possible to reproduce the morphology and kinematics of these encounters and that a good match 
strongly constrains the encounter parameters. 

Subject headings: galaxies: interactions - galaxies: kinematics and dynamics - methods: N-body simulations 



1. INTRODUCTION 

The diverse morphological and kinematic features of in- 
Q teracting disk galaxies have a simple dynamical explana- 
$_i tion: galactic bridges and tails (Toomre & Toomre 1972), 
^ rings (Lynds & Toomre 1976; Theys & Spiegel 1977) and re- 
Ci lated structures result when ordinary galactic disks experience 
~~' strong tides in close encounters. Tides also cause interact- 
, ing galaxies to merge by inexorably transferring energy and 
momentum from relative motion to internal degrees of free- 
dom (Toomre 1977; White 1978; Barnes 1988). With such a 
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straightforward physical basis, one might expect that dynami- 
cal modeling of interacting galaxies would be relatively easy. 
However, it's very time-consuming to explore the large pa- 
rameter space required to describe a galaxy collision and find 
a good match to the kinematics and morphology of a specific 
system. In addition, it's never been entirely clear that a good 
match yields a unique or physically meaningful model. 

In this paper we develop and test an efficient methodology 
to model the observable morphology and kinematics of pairs 
of interacting disk galaxies. Empirical tests show that the re- 
sulting models can be used to make strong inferences about 
the systems they match. While we do not consider minor 
mergers and interactions in this paper, our methodology can 
easily be extended to treat such encounters. 

At first glance, it seems all too easy to model interacting 
galaxies - and impossible to do so with any degree of confi- 
dence. The dynamical state of a galactic collision is described 
by a phase-space distribution function, /(r,v), which gives 
the mass density at position r and velocity v. In contrast, ob- 
servations of a specific component c (e.g., neutral hydrogen) 
yield a data cube Fc{X ,Y,V), which represents the distribu- 
tion of that component at each point {X,Y) on the plane of 
the sky as a function of line-of-sight velocity V. Since /(r , v) 
depends on six variables, while Fc{X,Y,V) depends on only 
three, it appears that observations do not provide enough in- 

'Xhis software is available at http://www.ifa.hawaii.edu/faculty/barnes/ 
research/identikit/ . 



formation. Put simply, there are an infinite number of differ- 
ent 6-D distribution functions consistent with any given 3-D 
data cube. 

On further reflection the problem is not quite as hopeless 
as it appears. This is because a typical galaxy merger begin 
with a tidal encounter between two normal, fairly symmetric 
spirals. Galaxies are scrambled as they merge, but the stars 
and dark matter which constitute most of their mass evolve 
collisionlessly. The fundamental dynamical equation. 



df df df ^ 

— + V • — ^ = 

dt dr dr d\ ' 



(1) 



where <I> is the gravitational potential, is fully reversible (e.g. 
van Albada & van Gorkom 1977); thus in some sense the orig- 
inal galaxies are still there, imposing a hidden symmetry on 
the dynamical state of a merging system. So we need not 
consider all possible distribution functions consistent with a 
given data cube; only a very small subset of these functions 
can possibly result from an encounter between two normal 
disk galaxies. 

In practice, mergers are modeled by guessing initial condi- 
tions, numerically simulating the ensuing collision, and com- 
paring the result to the morphology and kinematics of the sys- 
tem one wants to model. If the model fails to match the ob- 
servations, go back and guess again until the results are satis- 
factory. Most of the guess-work focuses on selecting the disk 
orientations, typically specified by angles (/i , (Oi ) and (/2, (Oi) 
for disks 1 and 2 (Toomre & Toomre 1972); also needed are 
the eccentricity e and pericentric separation p of the initial or- 
bit, as well as the galactic mass ratio ji. So seven parameters 
are needed to specify the initial conditions for an encounter of 
two axisymmetric disk galaxies - not counting the parameters 
used to specify their internal structures ' . 

' Selecting the correct internal structures is a separate problem, and one 
largely beyond the scope of this paper. Under some fairly general assump- 
tions, the internal structure of an axisymmetric galaxy may be described by 
a distribution function / = f(E,J-i^J^) depending on the energy E, angular 
momentum about the symmetry axis y^, and a third integral of motion 73. 
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Once a simulation has been run, one must select another 
nine parameters when comparing the results to observational 
data: a time t since pericenter, a viewing direction given, 
for example, by angles (0x,0y,0z), scale factors ^ and Y 
for length and velocity^, and a center-of-mass position on 
the plane of the sky (Xcm,Jcm) and velocity Vcm- All told, 
a minimum of sixteen parameters are needed to completely 
specify the initial conditions, time, point of view, and scale 
of a merger model. This plethora of parameters has long 
posed a challenge for systematic surveys of galactic colUsions 
(Toomre & Toomre 1972; Farouki & Shapiro 1982; Wallin & 
Stuart 1992; Howard et al. 1993; Barnes 1998; Naab & Burk- 
ert 2003). The problem we address here is slightly differ- 
ent - instead of trying to survey the entire parameter space, 
we want to navigate toward a solution matching the morphol- 
ogy and kinematics of a given interacting system. Intuition 
and prior experience can guide this process by narrowing the 
range of parameter space explored. Nonetheless, given the 
size and complexity of this parameter space, it's not surpris- 
ing that many simulations must be run to attempt a match, 
or that detailed models of galactic colUsions are not easy to 
produce. 



2. "IDENTIKIT" METHODOLOGY 

While self-consistent simulations are useful to finalize dy- 
namical models of tidally interacting galaxies, here we simu- 
late galactic disks with test particles. Test particles have a long 
history (Pfleiderer & Seidentopf 1961; Toomre & Toomre 
1972; Combes et al. 1980; Hernquist & Quinn 1987; Wallin 
& Stuart 1992) and nicely reproduce features such as bridges, 
tails, and shells which develop with little direct influence 
from self-gravity. To include orbital decay, which is crucial 
in modeling the more advanced stages of galaxy encounters 
and mergers, the test particles may be used to estimate the 
drag on the central masses (e.g., Toomre & Toomre 1972; 
Borne 1984; Quinn & Goodman 1986). However, orbit decay 
is largely driven by tidal interactions of galaxy halos (e.g., 
Toomre 1977; White 1978; Barnes 1988), and it's relatively 
easy to compute the self-consistent interaction of two halos 
using N-body simulations. Our initial approach was therefore 
to represent the mass of each galaxy with a spherical distribu- 
tion of massive particles; in each of these spheres, we embed- 
ded multiple disks of test particles, and decided which disk to 
display after running the simulations^. 

Building on this idea, and taking advantage of the faster 
processors now available, we have replaced these discrete 
collections of disks with spherical swarms of test particles 
moving on circular orbits - in effect, populating each galaxy 
model with all possible disks. 

The modeling procedure begins by selecting mass models 
for the two galaxies, thereby fixing the mass ratio ii. Our 
models include a bulge, a disk, and a halo; these compo- 
nents have cumulative mass profiles mb(r), mi{r), and mh(r), 
respectively, where mc{r) is the total mass in component c 
within radius r. For each galaxy, we compute the total mass 

Formally speaking, an infinite number of parameters are needed to specify 

such a function. 

' If tlie simulation is conducted in physical units then these parameters are 
not necessary - but additional parameters are required to describe the initial 
conditions, so the total parameter count is unchanged. 

^ In forensic investigations, "Identikit" is one of several systems used to 
construct portraits by selecting from a menu of facial features. Our approach 
is analogous. 



profile 

m{r)=mb{r)+mi{r)+nn,{r), (2) 

calculate the corresponding isotropic distribution function us- 
ing Eddington's (1916) formula (e.g., Binney & Tremaine 
1987, p. 236), and construct a spherical N-body realization 
of this profile using A^sphr equal-mass particles. Details of this 
procedure are given in Appendix A. 

Each of these spherical N-body realizations is then loaded 
with Attest test particles moving on circular orbits. We pick the 
orientation of the orbit of test particle / by randomly drawing 
its normalized angular momentum s, from a uniform distribu- 
tion on the unit sphere S'^. The radial distribution of the test 
particles may be chosen at will. The simplest choice is to use 
a radial distribution following the cumulative profile of the 
disk, mi{r); with this approach, the test particles representing 
a disk with a normaUzed spin vector Sd are those with 

1-Sd-S,<(7, (3) 

where (7 <C 1 is a tolerance parameter proportional to the num- 
ber of particles selected. In practice, however, this places 
many disk particles at small radii where they are largely 
immune to tides. The sampUng at large radii can be im- 
proved by radially biasing the test particle distribution. We 
do this by multiplying the local disk particle density, PA{r) = 
(Anr^)^^dmi/dr, by a factor of r^, and replacing (3) with 

1 -Sd-S; < a/max(^,,rmin)^. (4) 

Here qt is the initial orbital radius of particle ;, and r^i^ is a 
parameter which keeps (4) from diverging for small qt. Par- 
ticles selected using (4) follow the original disk distribution 
down to radius fmin, at smaller radii the disk is undersampled, 
but this has little effect if rmin is small. 

Finally, two such configurations are placed on a relative or- 
bit with a given pericentric separation p and eccentricity e and 
followed until they merge; we save particle positions and ve- 
locities every few time-steps, creating a data-base of several 
hundred frames tracing the system's history from start to fin- 
ish. This data can then be used to approximate any encounter 
with the chosen /i, p, and e. 

Identikit software includes interactive routines allowing the 
user to select the disk orientations, viewing direction, scale 
factors, and centroid positions; the resulting test-particle co- 
ordinates are instantly projected on the (X,y), {X,V), {V,Y), 
and (X,Z) planes. The user can also step forward or back- 
ward in time, switch between data-bases created using differ- 
ent values for p, jj., or e, and vary the tolerance parameter a. 
In addition, an observational data cube F{X,Y,V) for a spe- 
cific system to be matched may be projected onto the {X,Y), 
{X,V), and {V,Y) planes, typically using grey-scale or con- 
tour images. A lap-top computer can easily store the neces- 
sary data and supply the modest processing power required 
when searching for a match. Fig. 1 presents an example; the 
encounter and viewing parameters have been adjusted so the 
test particles (points) closely match the data cube (grey-scale). 

3. EXPERIMENTAL VALIDATION 

The Identikit system is designed with the primary goal of 
analyzing observations. However, it's not straightforward to 
perform empirical tests using observational data. For one 
thing, much of the available data is rather heterogeneous; 
data-sets with uniformly high resolution and signal-to-noise 
ratio are not easy to obtain. For another, we don't know the 
true initial conditions and viewing parameters for more than a 



Fig. 1. — Identikit match to a pair of merging galaxies. Grey-scale images represent the data cube F{X,Y,V) of the system, while points represent the 
Identikit model. Top-left, top-right, and bottom-left quadrants show the data cube and the model projected on the (X,Y), (V,Y), and {X,V) planes, respectively; 
bottom-right quadrant shows the model projected on the {X,Z) plane. 



handful of the mergers that have been observed, so we don't 
have any simple way to validate the results of our modeling. 
To determine if Identikit can actually reconstruct galactic en- 
counters from the information contained in data-cubes, we 
tested it on an artificial data set of 36 self-consistent disk 
galaxy merger simulations with random orientations, times 
since first passage, viewing directions, and scale factors. 

3.1. Artificial Merger Data 

Our disk galaxy model has a spherical bulge (Hemquist 
1990a) containing 5% of the mass, an exponential/isothermal 
disk (de Vaucouleurs 1959a,b; Freeman 1970; van der Kruit 
& Searle 1981) containing another 15%, and a spherical dark 
halo (Navarro et al. 1996) containing the remaining 80%. The 
density profiles for these components are 

Pb(r)ocr-l(r-fab)"\ 

Pd{q,z) « e-^/^o sech2(z/zd) , (5) 

Ph(r)ocr-i(r-Fah)"^, 

where is the scale length of the bulge, q = \/x^ is 
the cylindrical radius, ad is the scale length of the disk, zd is 



the scale height of the disk, and is the scale length of the 
halo. Each galaxy was realized using a total of = 131072 
particles. The simulations used natural units with Newton's 
constant G = 1. In these units, the galaxy model has total 
mass m = 1 .25 and half-mass radius rn,cd — 0.532. The disk's 
scale length Od — 1/12, and the median circular velocity of 
the disk material is Vmed — 1-23; at a radius of 3ad the orbital 
period is ?orb — 1-23. 

We restricted our artificial data set to equal-mass (n = I) 
encounters with parabolic initial orbits (e = 1); for each orbit, 
the pericentric separation p was drawn from a uniform distri- 
bution in the range [0.05,0.5] = [0.6,6]fld- The instant when 
this idealized two-body orbit reaches pericenter defines t = 0; 
times ? < are before pericenter, while times ? > are after 
pericenter. We adopt a coordinate system in which the orbital 
angular momentum vector is parallel to the z axis. The nor- 
malized spin vector § of each disk was chosen from a uniform 
distribution on the unit sphere S^; in practice, the inclination i 
was chosen by drawing cos(/) = § • z from a uniform distribu- 
tion in the range [—1,1], and the argument (0 was chosen from 
a uniform distribution in the range [0°, 360°]. Further details 



4 



on the galaxy models and the merger simulations are given in 
Appendix B. 

For each of the 36 merger simulations we chose a random 
time between first and second pericenter, rescaled the system 
by random factors in length and velocity, and "observed" it 
from a random direction. We first determined the relative orbit 
of each pair of galaxies, using the most tightly-bound 2048 
particles in each bulge to measure galactic positions. Let t[ ~ 
and t2 be times of first and second pericenter, respectively; 
the random time t was drawn from a uniform distribution in 
the range [?i,?2] and rounded down to the nearest available 
output time. 

Next, we selected scale factors ^ and Y for length and 
velocity, respectively. These were chosen so that the galaxy 
models obey a mass-radius- velocity relation of the form M oc 
oc y4 with a small amount of scatter (c.f. Tully & Fisher 
1977). Let ^ be drawn from a uniform distribution in the 
range [—0.5,0.5], and and g2 be drawn from a Gaussian 
distribution with zero mean and unit dispersion; then 

^ = 10«/2 100 05«l , -r = 10^/4 100 05^2 . (6) 

Finally, we chose a random viewing direction Z from a uni- 
form distribution on the unit sphere S'^. We drew a second 
vector Xo from the same distribution, and set X = Xq — Z(Xo • 
Z) and Y = Z X X. These vectors and scale factors were used 
to map the position r, and velocity v, of each particle i to data- 
cube coordinates: 

Xi = ^X-ri, Yi = ££\-Xi, V; = rZ-Vi. (7) 

Particles from the disks of the two galaxies, transformed 
to (X,,y,, V,) coordinates, provide an N-body representation 
of a data cube Fi\{X^Y,Y) for the disk material. Such data 
is roughly comparable to the neutral hydrogen data-cubes 
FHi(o!,5,yios) available for many interacting galaxies (e.g., 
Hibbard et al. 2001). The simulated data has better reso- 
lution than most observational data-sets and is free of noise 
and interferometric artifacts; moreover, our simulations used 
collisionless particles instead of neutial gas. We could have 
run random mergers with gas to improve the correspondence 
between the simulations and real observational data, but the 
computing time required for a large suite of simulations with 
gas is non-trivial. Fortunately, collisionless simulations do a 
good job of reproducing the tidal features commonly detected 
in HI since the latter usually evolve balhstically once tidally 
extracted from their parent galaxies. 

To present the simulated data in the form required for Iden- 
tikit matching, we projected the disk particles for each of our 
36 mergers on the (X,y), and (VjF) planes; gridded 

particle distributions were Ughtiy smoothed to produce grey- 
scale images. We also used tightly-bound particles from the 
bulge of each galaxy to determine its position and line-of- 
sight velocity; the results were plotted on top of the grey-scale 
images. Our images are thus analogous to HI maps supple- 
mented with accurate nuclear coordinates and velocities. 

The entire procedure outUned above, including both the 
generation of the simulations and the selection of viewing pa- 
rameters, was performed by automated scripts without human 
intervention; we did not know the actual values of any pa- 
rameters except e and \i. The resulting sample of merging 
galaxies, shown in Fig. 2, possess a variety of morphologies; 
only a subset display the "double tails" characteristic of the 
best-known mergers (Toomre 1977). 



3.2. Identikit Matching 

We prepared a series of eight Identikit simulations spanning 
a range of pericentric separations. Each simulation contained 
two identical (/i = 1) configurations of Attest = 262144 test 
particles and A^sphr = 81920 massive particles; the mass model 
used to set up the massive particle distribution was a spheri- 
cal version of the one used in the random mergers. These 
configurations were placed on parabolic (e = 1) relative or- 
bits with pericentric separations p= 1/16, 2/16, 3/16, 4/16, 
5/16, 6/16, 7/16, and 8/16, starting at f = —2 time units be- 
fore first pericenter, and followed until ? = 8, by which time 
even the widest passage had merged. 

We used these simulations and the Identikit software to 
fit each of the random mergers by interactively matching 
the "observed" (X,y), (X,Y), and (y,i') projections with 
test particles. The modeling process usually began with 
rough guesses for the viewing direction, time since pericen- 
ter, and pericentric separation. A variety of clues guided these 
guesses. For example, short but pronounced tidal features 
point to a recent tidal encounter, while long but attenuated 
features suggest an older passage; in later stages, loops asso- 
ciated with tails show that material has started falling back. 
Likewise, if the galaxies display a large separation in pro- 
jected velocity then the sight-Une must be close to the orbital 
plane; conversely, a small difference in projected velocity im- 
plies either that the system is observed near apocenter or that 
the relative velocity vector is roughly perpendicular to the line 
of sight. Finally, other things being equal, closer passages 
generally yield stronger and more dramatic tidal features. 

The next step was to adjust the orientations of the two disks, 
attempting to roughly match the morphology and kinematics 
of the system. This generally suggested further modifications 
to the viewing direction, separation, and time, as well as the 
scale factors and center-of-mass position and velocity. Fur- 
ther adjustment of all parameters continued until a satisfac- 
tory match was obtained or exhaustion set in. Our criteria for 
a satisfactory match were somewhat subjective^; we placed a 
good deal of weight on matching tidal features (e.g.. Fig. 1), 
while recognizing that test particles can't accurately repro- 
duce structures - such as tidally-induced spirals - which de- 
pend on self-gravity. Each match typically took a few hours, 
and the entire set of 36 random mergers was matched in about 
one month; for comparison, a match to the NGC 7252 merger 
remnant (Hibbard & Mihos 1995) took 74 N-body runs over 
a three month time period, while a match to the NGC 4676 
system (Barnes 2004) took ^ 30 runs over two months. 

After changing the viewing direction, time since pericen- 
ter or orbital parameters, it's usually necessary to reposition 
the centers of the models on top of the actual positions by 
adjusting the rotation about the viewing axis dx, scale fac- 
tor 5£ , and center of mass position {^can^cm)- We therefore 
implemented an option to "lock" the centers; when this op- 
tion is invoked, Bx, and (Xcm,Jcm) are recalculated on the 
fly, keeping the projected positions of the models invariant as 
other parameters are changed. Locking works quite well when 
the two galaxies are well-separated on the plane; it's 

less useful, and can be downright counter-productive, when 
the centers appear close together In fine-tuning a nearly fi- 
nal match we sometimes found it useful to unlock the centers, 
trading off slight misahgnments in central positions for im- 

It's not trivial to evaluate matches quantitatively; for more on this, see 
§4.2.1. 
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Fig. 2. — Sky-plane (X,y) projections of the 36 merging encounters used to test the Identikit procedure. Note the range of sizes and morphologies. The 
number of each system appears to its upper left; color indicates Identikit fit quality (see § 3.2), with good fits in red, fair fits in black, and poor fits in blue. 
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proved matches to tidal features. 

Fig. 1, which shows our match to object 23 (see Fig. 2), 
illustrates many aspects of the matching process. From the 
start, it seemed likely that the viewing direction would be 
fairly close to the orbital plane, since the two galaxies have 
rather different systemic velocities. The galaxy on the lower 
right of the {X,Y) projection appears nearly edge-on, as in- 
dicated by its morphology and its rather large velocity range. 
Since its tidal features lie more or less in the same plane as the 
disk itself, it seemed plausible that this galaxy has a relatively 
small incUnation ii to the orbital plane, while its companion 
clearly has a higher inclination i2 and appears more face-on 
from our viewpoint. The dual-valued velocities along the tail 
of the edge-on disk, which produce the "hook"-shaped feature 
seen in the {X,V) and {V,Y) projections, suggested that this 
tail is actually quite extended, but viewed so as to double back 
on itself. 

For an initial match to this system, we tried a "middle- 
of-the-road" pericentric separation (p — 0.25); at a relatively 
early time {t = 0.56) we could roughly match the velocity dif- 
ference and some aspects of the morphology and kinematics, 
including the spiral morphology of the face-on disk and the 
hooked tail in the {V,Y) projection. However, other features 
of this initial match were less satisfactory. In the {X,V) pro- 
jection, the tail doubled back too soon, while in the {V,Y) 
projection, the bridge did not span the velocity range between 
the galaxies, falhng to the left of its ideal position. Moreover, 
the more face-on disk, while nicely rendered in the {X,Y) 
projection, populated regions of phase space which the {V,Y ) 
projection showed to be empty. Trial and error revealed that 
wider passages and later times could repair most of these de- 
fects; the solution shown in Fig. 1 uses a p = 0.5 passage 
viewed at f = 1 . The tail in this match, while more extended 
than it was initially, is still a bit too short. Times f > 1 yield 
longer tails, but the velocity difference between the galaxies 
becomes too small, and bridge particles falling through the 
more face-on disk over-populate a relatively sparse region of 
phase-space. The adopted solution is therefore a compromise 
between several competing factors. 

After comparing our Identikit models to the morphology 
and kinematics of all 36 random mergers, we subjectively 
graded the solutions as "good" (18 cases), "fair" (12 cases), 
or "poor" (6 cases); these grades are indicated in Fig. 2. Good 
matches, like the one in Fig. 1, strongly constrain the param- 
eters. Fair matches generally appear plausible but allow more 
latitude in selecting parameter values; this group included 
several systems with twin edge-on tidal tails. Poor matches 
could be divided into two groups: systems with weak and 
diffuse tidal features, typically resulting from very wide en- 
counters involving retrograde or highly inclined galaxies (ob- 
jects 16 and 32), and systems with pronounced but confusing 
tidal features (objects 1,6, 11, and 24). 

3.3. Results: Parameters 

With the Identikit solutions in hand, we compared their pa- 
rameter values to the true values used to generate the artifi- 
cial merger data. Fig. 3 uses an abstract cylindrical coordi- 
nate system to portray the parameter space explored in these 
solutions. The radial coordinate represents the initial orbit 
of the two galaxies; the fits discussed here parameterize the 
orbit by the pericentric separation p, since the eccentricity 
e and mass ratio n were fixed beforehand. The azimuthal 
coordinate represents the four angles (;i,a)i) and (/2,ti)2) re- 
quired to specify the initial orientations of the two disks. To- 
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Fig. 3. — An abstract representation of the parameter space of galaxy in- 
teractions. The radial coordinate represents the initial orbit, the azimuthal 
coordinate represents the disk orientations, and the vertical coordinate repre- 
sents the parameters chosen after a simulation is run. A conventional N-body 
simulation explores the parameter subspace represented by the dotted Une, 
while a single Identikit simulation can explore entire cyUndrical surface. 



gether, the radial and azimuthal coordinates of this abstract 
space completely specify the initial conditions for a galaxy 
interaction. The vertical coordinate represents the parameters 
selected after running a simulation: time since pericenter /, 
viewing angles (0x,^,fe), and scale factors (^,1^); here 
the center-of-mass parameters are omitted since their values 
are not discussed below. A conventional N-body simulation 
starts at a point on the horizontal plane and explores the pa- 
rameter subspace represented by the dotted line in this figure; 
a single Identikit simulation, in contrast, allows access to an 
entire cylindrical surface. 

Fig. 4 presents a scatter plot comparing inclinations /fit de- 
rived from the Identikit fits against the inclinations ;t,uc used 
in the random merger sample. Here and in subsequent plots, 
color indicates the grade of each model; note that the disks are 
not graded individually, so both disks in a given model receive 
the same grade even if one fits better than the other. The good 
fits (shown in red) fall quite close to the diagonal line across 
the entire range of the plot. The fair fits (black) display more 
scatter but track the same relationship. In contrast, the poor 
fits (blue) have a very different distribution: almost all have 
incUnations Jtme > 90°, and most fall quite far from the diago- 
nal. This plot supports a couple of useful inferences. First, our 
subjective grades, based on the overall appearance of the Iden- 
tikit models, correlate with |ifit — itme|; in other words, these 
grades mean something. Second, encounters involving disks 
with incUnations i > 90° are more difficult to model, presum- 
ably because the tidal features such encounters produce are 
less distinct and more ambiguous. Nonetheless, many high- 
inclination encounters were successfully modeled; other fac- 
tors evidently influence the outcome of the modeling process. 

Likewise, Fig. 5 compares fit and true values of the argu- 
ment to pericenter, CDfit and COtme- Since (0 becomes indeter- 
minate for inclinations near ; = 0° or / = 180°, we plot disks 
with 30° < 'true l£ 150° as filled circles, and disks outside this 
range as crosses. This plot shows that good fits yield argu- 
ments quite close to the true values, and the fair fits do nearly 
as well. At first sight it may seem that a few disks, repre- 
sented by the one good and three fair points in the upper left 
and lower right of the plot, yield discrepant values of Ofit, but 
this is an artifact of topology; 0) is a periodic coordinate, and 
when opposing edges of the plot are identified, these appar- 
ent outliers are not so far from cOgt = (Otrue- The poor fits, in 
contrast, genuinely scatter throughout the plot. Most of the 
good and fair fits with ixme < 30° or 150° < itme still yield 
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Fig. 4. — Estimates of disk inclination, r. The color of each data point 
shows quality of fit: red is good, black is fair, and blue is poor. The heavy 
dotted line represents perfect agreement ((jit = i'lrue); the light lines show ifn = 
W± 10°. 

reasonable values for (O, but these are a bit more scattered. 
For example, only one of the six good fits with /true 

in this 

range yields a coat within 10° of Otme; in contrast, just over 
half of the good fits with inclinations 30° < itme 

< 150° have 

|COfit-fiJtrue| < 10°. 

To examine our overall accuracy in determining initial disk 
orientations, we computed misalignments A^^m between true 
and fitted spin vectors for all 72 disks. Fig. 6 shows the cumu- 
lative distribution function of cos(Aspin) = s^t • Stme; as in pre- 
vious plots, color indicates the subjective grade of each fit. A 
perfect match yields cos(Aspm) — 1, while if Sm was uncorre- 
cted with Stme then cos(Aspin) would be uniformly distributed 
in the range [—1,1] and the plotted points would fall along 
a diagonal from lower left to upper right. The actual distri- 
bution is strongly peaked near cos(Aspin) ~ 1, with the good 
fits showing the smallest misalignments, the poor fits showing 
the largest misalignments, and the fair fits falling in between. 
This is consistent with the previous figures, since / and co are 
just angular coordinates for s. For the entire sample, the me- 
dian value is Aspin ~ 18°, while for the 36 disks in good fits 
the median is Aspin = 12°. 

In a similar fashion. Fig. 7 shows the distribution of the 
misalignment in viewing direction, Ayiew for all 36 fits. 
Here cos(Aview) = Zfjt • Ztme, perfect agreement again yields 
cos(Aview) = 1, and perfect ignorance would distribute points 
along a diagonal from lower left to upper right. Viewing di- 
rection is quite well determined; for the entire sample the me- 
dian Aview = 13°, while for the 18 good fits alone the median 
is only slightly smaller, Ayjew = 12°. Note that an error circle 
with a radius of 12° covers roughly 1% of the solid angle of a 
sphere; these fits are much better than educated guesses! 

Fig. 8 compares the actual time ftme since pericenter against 
the time fflt obtained from the Identikit fit. The scatter plot on 
the left shows that the fitted and true values are in good agree- 
ment, closely tracking each other throughout the entire range 
of times; there is no evidence of bias or systematic error, and 
the residuals appear to be random. The plot on the right shows 
the cumulative distribution of the fit/true ratio, ffit/ftme- The 



Fig. 5. — Estimates of disk argument, (O. Colors and dotted lines as in 
Fig. 4; crosses represent disks with | ('ime — 90° | > 60° . Note that opposing 
edges of this plot should be identified. 

symmetric appearance of this curve provides further evidence 
that f is accurately estimated by the Identikit models. Sub- 
jective fit quality appears to correlate with ffit/ftrue; of the ten 
points at the two extremes of the distribution, only one comes 
from a good fit. 

Fig. 9 compares fitted and true values of the pericentric sep- 
aration. Here the range of p values is rather small, and the 
fact that pfit can take on only eight discrete values is evident. 
There is a fair correlation between pfit and ptme, although the 
points show considerable scatter. The cumulative distribution 
of pnt/ptme plotted on the right is nonetheless fairly symmet- 
ric, and the median value of pnt/ptmc is very close to unity. 
Curiously, there's not much sign that the grades assigned the 
models correlate with pf^i/ ptme; fair and poor fits appear inter- 
spersed with good ones throughout most of the distribution. 

Identikit estimates of the length and velocity scale factors, 
^ and y, are compared with their true values in Figs. 10 
and 11, respectively. The length scale ^ is well-determined; 
the plot on the left shows points scattered about the diago- 
nal line, while the cumulative distribution on the right shows 
a narrow spread with a median value of ^ntf-^tme close to 
unity. In contrast. Identikit estimates of the velocity scale 
Y show a small but very definite bias; the distribution of 
"/fit/l^rue is narrow, but clearly offset from unity. Tests with 
isolated disks suggest that the absence of random motion 
in the Identikit test-particle disks is responsible; see Ap- 
pendix C. 

Table 1 lists statistics for the ratios of fitted to true values 
for f , p, and T^; the physical parameters T and P will be 
discussed in the next section. Fit/true ratios for t and ^ have 
fairly narrow distributions centered on unity, confirming that 
these parameters are determined accurately and without bias. 
The fit/true distribution for p, while somewhat broader, is also 
centered on unity. The fit/true distribution for the velocity 
scale factor Y is quite narrow, but the median value is ~ 10% 
too high, showing again that a small bias is present in fitting 

r. 
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Fig. 6. — Cumulative distribution function for cos(Aspin). Color shows fit 
quality. 



Fig. 7. — Cumulative distribution function for cos(Avicw)- Color shows fit 
quality. 
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Fig. 8. — Estimates of the dimensionless time since pericenter, t. Left: scatter plot of log^jt against logJtm 
log?flt/?tme. In both plots, color shows fit quality, and the dotted line represents perfect agreement {tfu = ftme)- 



Right: cumulative distribution function of 



3.4. Results: Residuals 

The errors determined by comparing the Identikit mod- 
els with the actual mergers define six independent residu- 
als: viewing dkection (Ayiew), spin direction (Aspin), time 
since pericenter (ffit/ftme). separation at pericenter {pat/ Ptme). 
length scale (-Sffit/^tme), and velocity scale (1^t/Xrue)- If 
each of these residuals is plotted against the others, a total 
of fifteen (6x5/2) potential relationships can be examined. 
Most of these plots show no measurable correlation; only the 
four in Fig. 12 are significant (correlation coefficient > 0.5). 

As shown in the upper left panel of Fig. 12, errors in disk 
orientation Aspin are correlated with errors in viewing direc- 
tion Ayiew- This correlation, which is largely driven by the 
poor fits, is not hard to explain. Once a viewing direction has 



been selected, the next step is usually to adjust the disk ori- 
entations; if clear guidance from tidal features is lacking, the 
best one can do is to match each disk's position angle and ap- 
parent tilt with respect to the line of sight. A poor choice for 
the viewing direction will induce comparable errors in spin 
direction, as seen here. 

The plot in the upper right panel shows that residuals in the 
time since pericenter t correlate with residuals in pericentric 
separation p. This correlation also has a simple explanation; 
wider passages evolve and merge more slowly, so if ffit is for 
some reason overestimated during the matching process then 
a larger pf^t can partly compensate for this error. As an ex- 
treme example, suppose tfn was set so high that an encounter 
with the correct pf^ would already have merged by this time; 
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Fig. 9. — Estimates of the dimensionless pericentric separation, p. Left: scatter plot of pflt against piruc- Right: cumulative distribution function of Pflt/pime- 
Colors and dotted lines as in Fig. 8. 
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Fig. 10. — Estimates of the length scale factor, Left: scatter plot of log^jt against log^tiuc- Right: cumulative distribution function of ^flt/^tiue- Colors 
and dotted lines as in Fig. 8. 



by selecting a larger value for pf^, the merger can be delayed 
until some time after ffit, improving the match to the morphol- 
ogy- 

A clear correlation between the residuals in t and of 
the form J^fn/^tme — (ffit/ftrue) ' is shown in the lower right 
panel of Fig. 12. Such a correlation arises naturally in cases 
where fits are strongly constrained by well-developed tidal 
tails. In a proper tail, material at the tip is moving fast enough 
to escape; at late times this material will have asymptotically 
constant velocity (Toomre & Toomre 1972), implying that a 
tail's length grows in proportion to its age (Schweizer 1977). 
If ffit is overestimated, tidal tails will be longer in direct pro- 
portion, and a smaller value of will be needed to rescale 



them to their correct physical lengths. It's likely that similar 
considerations also apply when other kinds of tidal features 
provide the primary constraints; this is consistent with the re- 
sults presented here, which show that the good and fair fits 
fall close to ^fn/ J^ime — (ffit/ftme) while the poor ones are 
more scattered. 

Finally, we conjecture that the correlation between residu- 
als of p and ^ shown in the lower left panel is induced by the 
ffit/ftrue-Pfit/ptme and ffit/ftrue--%t/^trae coiTclations already 
described. Of the three, the ;?fit/;?true--%t/-S?true correlation 
shows the most scatter, and its outliers tend to match those in 
the ffit/ftrue-Pfit//'true Correlation. 

An interesting consequence of these correlations emerges 
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Fig. 1 1 . — Estimates of the velocity scale factor, "V . Left: scatter plot of log against logXrue- Right: cumulative distribution function of ^t/Xme- Colors 
and dotted lines as in Fig. 8. 
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Fig. 12. — Scatter-plots of statistically significant correlations between 
residuals in Identikit models. Color indicates quality of fit. The plot on the 
upper left contains twice as many points as the others because Ajpi^ is plotted 
for each disk. 



when the dimensionless model parameters t and p are com- 
bined with the length and velocity scale factors ^ and . 
In matching an Identikit model to observational data, X and 
Y might be given in units of kpc, and V might be given in 
units of kms^'; the scale factors ^ and 'V would then have 
units of kpc and kms^^ respectively, and define a transfor- 
mation from dimensionless model data to real physical values. 
The physical pericentric separation is P = ^p, while the time 
since pericenter is T = /y)t. A glance at the two lower 
plots in Fig. 12 suggests that these correlations may actually 
reduce errors in estimates of the physical parameters T and P. 



Table 1 supports this; the uncertainty in P, as indicated by the 
ratio of to quartile values, is about half the uncertainty 
in p. For T the improvement is not as striking; the ratio of 
^rd 2 St quartiles is slightly smaller than the ratio for f , but it 
appears that ^ introduces some additional scatter and a small 
bias. 

3.5. Poor Matches 

For balance, it's worth taking a closer look at some of the 
less successful matches. We judged six out of 36 of our solu- 
tions to be poor fits to the "observed" data. As already noted, 
two of the systems with poor matches had very diffuse tidal 
features, while four had fairly strong features. Our difficulties 
in modeling the latter are somewhat surprising; why, given 
the information which must be present in strong tidal features, 
didn't the models turn out better? 

Fig. 13 presents three versions of one of the more difficult 
systems (object 6 in Fig. 2). For ease of comparison, only 
the {X,Y) and [X^V) planes are shown; the (y,F) views are 
harder to interpret since the galaxies partly overlap. The left- 
hand version shows an Identikit model with the correct disk 
orientations, viewing direction, and velocity scale "dialed in" 
by hand. Only a finite set of pericentric separations are avail- 
able in Identikit, so it was not possible to exactly match the 
pericentric separation /Jtme = 0.223; on the whole, pfn = 3/16 
gave a better match than pflt = 1/4 although the latter is nu- 
merically closer This closer approach favored a somewhat 
different length scale (^fit = 0.644 vs. ^true = 0.590) and a 
slightly earlier time (/fit = 1.62 vs. ftme = 1-69); both of these 
adjustments are qualitatively consistent with the residual cor- 
relations discussed in the last section. 

The "correct" solution on the left in Fig. 13 does a poor job 
of matching the morphology (top). It does better at matching 
the kinematics (bottom), although some fairly large regions of 
phase space are quite under-populated, and one of the spiral 
arms in the more face-on disk produces a badly-placed streak 
of particles across the {X,V) projection. By our subjective 
criteria, this solution is a poor fit. Since the disk orientations 





Fig. 13. — Three models in search of a match. Grey-scale images show the "observed" data (object 6 in Fig. 2), while points represent models. Top and bottom 
rows show projections on {X,Y) and {X,V) planes, respectively. Left: the "correct" Identikit solution. Middle: a match to the shape and velocities of the more 
edge-on disk. Right: a match to the spiral features of the more face-on disk. 



and viewing direction exactly match those used to generate 
the grey-scale images, this mismatch must be largely due to 
the lack of disk self-gravity in the Identikit simulations. In 
particular, the spiral morphology of the more face-on disk is 
very poorly reproduced; with hindsight, we recognize that this 
spiral probably owes its form to self-gravity. 

The middle and right parts of Fig. 13 present two of the 
solutions we found while trying to match this system. Es- 
sentially, each matches one galaxy while failing to match the 
other. The middle solution does a plausible job of matching 
the morphology and kinematics of the more edge-on disk, us- 
ing a close passage {pfn = 1/ 16) at an early time (fgt = 1-12). 
Although it does a very poor job with the more face-on disk, 
this solution is actually not too bad - the viewing angle is off 
by 47°, but the disks are within 30° of their correct orien- 
tations. The right-hand solution approximates the morphol- 
ogy and kinematics of the more face-on disk, using a wider 
passage (pm = 3/16) at the same early time. Ironically, al- 
though the decent-looking fit to the face-on disk and some- 
what plausible match to the edge-on disk led us to adopt this 
solution as the best of a bad lot, it's actually a much less ac- 
curate fit than the middle one - the viewing angle is off by 
95°, while the face-on and edge-on disks are misaligned by 
109° and 94°, respectively! The splattering of particles be- 
tween the galaxies in the (X,V) projection, populating a re- 
gion which should have been completely empty, was the crit- 
ical flaw which saved us from assigning a "fair" grade to this 
solution. 

The other three systems which we failed to fit despite their 
strong tidal features are a mixed lot. In object 24, much as 



in the example just discussed, the "correct" Identikit solution 
doesn't match the morphology of either galaxy. In object 11, 
the correct solution matches one disk but fails to reproduce 
the morphology and kinematics of its partner Finally, ob- 
ject 1 is reproduced quite accurately by the correct Identikit 
model. Our failure with object 1 was largely due to insuffi- 
cient patience, while with object 24 the mismatch was so bad 
that the correct solution does not stand out when compared to 
other possible fits. Object 11 is intermediate; had we stumbled 
across the correct match to one disk we might have recognized 
the overall plausibility of this solution. In sum, human limi- 
tations and lack of self-gravity contribute about equally to the 
four poor matches of systems with strong tidal features. 

4. DISCUSSION 

Under somewhat idealized conditions, the Identikit 
methodology can recover the key parameters of a galactic en- 
counter from a single data-cube. If the method can also be 
applied to real data, it will be a powerful tool for interpreting 
observations and reconstructing the dynamical histories of in- 
teracting galaxies. But before going further, some limitations 
should be discussed. 

4.1. Limitations 

First, as noted above. Identikit simulations are not com- 
pletely self-consistent. Halos and bulges are treated self- 
consistently, but disks - of necessity - must be modeled with 
test particles. Disk structures requiring self-gravity, includ- 
ing bars and swing-amplified spirals, will not be reproduced. 
Consequently, these features can't be used to match models to 
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observations. Moreover, orbital decay of Identikit models is 
independent of disk orientation, whereas direct passages are 
expected to decay faster than retrograde ones (White 1979; 
Barnes 1992). For a pair of galaxies approaching their sec- 
ond passage, the estimated time tfn since first pericenter may 
be off by as much as ^ 20%; in applications requiring an ac- 
curate estimate of f^t. Identikit models should be followed up 
with self-consistent simulations. 

Second, while the test particles we used to model disks 
are collisionless, most kinematic tracers foUow a specific 
phase of the interstellar medium - for example, HI, Ha, or 
CO. Collisionless particles can approximate gas dynamics, 
but only if the gas moves ballistically; streams of particles 
freely interpenetrate, whereas gas will be deflected if it en- 
counters shocks. Examples include mass transfer via gen- 
uine bridges formed in low-inclination encounters (Toomre 
& Toomre 1972) and "splash bridges" due to hydrodynamic 
forces in interpenetrating encounters (Struck 1997). More 
subtle are the modest star-gas offsets seen in low-inclination 
encounters with extended gas disks (Mihos 2001); it appears 
that the gas offset from the stars does foUow ballistic trajec- 
tories, while some gas initially associated with the stars dis- 
sipates and falls back. Additionally, shocks may change the 
physical state of the gas; in particular, dense molecular gas 
is often associated with material which has undergone sig- 
nificant dissipation and hence cannot be modeled colUsion- 
lessly. But for most encounter geometries, gas in tidal fea- 
tures should be well-approximated by collisionless particles; 
exceptions can be recognized and allowed for in the fitting 
process. Stars, especially populations pre-dating the onset of 
an encounter, could be a useful complement to gas-phase trac- 
ers; absorption-line spectroscopy would require prohibitively 
large amounts of telescope time, but individual planetary neb- 
ulae are already providing kinematic data for nearby systems 
(e.g., Hui et al. 1995; Durrell et al. 2003). 

Third, two parameters were not included in the present ex- 
periment: the orbital eccentricity e and the mass ratio [i. In 
principle it's straightforward to include these parameters in 
the fitting process, although doing so increases the number of 
Identikit models needed. Both parameters have a priori con- 
straints - orbits with e ~ 1 are favored theoretically (Toomre 
& Toomre 1972; but see Khochfar & Burkert 2006), while n 
can sometimes be estimated photometrically assuming con- 
stant M/L ratios - so it seemed reasonable to exclude them 
from the initial experiments. Putting these parameters in play 
would probably increase the scatter in our fits; moreover, the 
three-way correlation between residuals of t, p, and ^ (§ 3.4) 
might expand to involve residuals of e and jti as well. Within 
the context of the models examined here, there's scant rea- 
son to expect much cross-talk between these parameters and 
others; in particular, disk orientations and viewing angles are 
robustly constrained by tidal features and should be well de- 
termined even if e and jU are included in the fits. 

Fourth, real galaxies have a range of rotation curve shapes, 
reflecting a diversity of mass profiles (e.g., Casertano & van 
Gorkom 1991; Catinella et al. 2006); a single mass model is 
too limited. The fact that we used the same mass model for 
our artificial data and for the Identikit simulations no doubt 
helped to reduce the uncertainties in our solutions. There's 
no reason why Identikit simulations cannot include a vari- 
ety of mass models, although this will increase the number 
of choices to be made in modeling a galactic collision. On the 
other hand, if such simulations can discriminate between dif- 
ferent mass models then they would provide a way to analyze 



the structure of disk galaxies. 

4.2. Previous Studies 

Models of interacting galaxies have a long history. Toomre 
& Toomre (1972) presented test-particle models of four sys- 
tems: Arp 295, M 51, NGC 4676, and NGC 4038/9; the lat- 
ter three have been revisited time and again by other work- 
ers. Kinematic information was initially scarce and of uneven 
quality, so early modeling attempts focused on reproducing 
the optical morphology. Better velocity data has given kine- 
matics a more substantial role in more recent modeling ef- 
forts. At the same time, faster computers and N-body algo- 
rithms have enabled researchers to construct models incorpo- 
rating self-gravity. 

The methodology adopted by Hibbard & Mihos (1995) in 
their model of the well-known merger remnant NGC 7252 in- 
cludes several key practices found in other successful models: 

1. Detafled velocity information, in the form of HI data 
(Hibbard et al. 1994), was available to constrain the 
model. 

2. Astrophysical arguments were used to estimate several 
critical parameters - specifically, the mass ratio jU, the 
initial orbital eccentricity e, and the pericentric separa- 
tion p. 

3. The model focused on reproducing the large-scale and 
morphology and kinematics of the tidal tails; these fea- 
tures evolved ballistically since first passage and there- 
fore carry a memory of the initial encounter. 

4. Test-particle models, with rigid galaxy potentials con- 
strained to follow realistic merger trajectories, were 
used to narrow down the range of parameter space. 

5 . Fully self-consistent N-body models were used to refine 
the final model. 

6. Simulation particles were plotted over orthogonal pro- 
jections of the data cube to show that the final model 
reproduced both the morphology and the kinematics of 
NGC 7252. 

Table 2 lists some interacting disk galaxies with dynami- 
cal models incorporating significant kinematical constraints. 
The progenitors of these systems span a range of mass ratios 
and morphological types. "Sh-S" systems involve two disk 
galaxies of roughly comparable mass, both generally display- 
ing significant tidal features. Most of these pairs are observed 
between first and second passage; NGC 7252 is the only 
completed merger. As Struck et al. (2005) note in modeling 
NGC 2207, earlier stages are generally easier to fit. In "S-nd" 
encounters a disk galaxy is perturbed by a smaller companion, 
while in "Eh-S" systems the disk is disturbed by an elliptical 
of comparable mass. Finally, in "ring" galaxies a compan- 
ion has plunged almost perpendicularly through a disk galaxy 
(Appleton & Struck-MarceU 1996); tiiese systems are rela- 
tively straightforward to model since their geometry is fairly 
simple. 

To varying degrees, the studies in Table 2 all followed the 
methodology used by Hibbard & Mihos (1995). Most had ac- 
cess to detailed velocity information, usually obtained by HI 
interferometry or Ha Fabry-Perot imaging, although a few 
models were based on long-sUt spectroscopy. Adopted orbital 
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eccentricities reflect a range of assumptions, not all equally 
plausible; for most systems, orbits with e ^ I seem more 
likely since e < 1 begs the question of what happened on the 
previous passage. Both test-particle and self-consistent tech- 
niques were used. Not all models were refined using fully 
self-consistent simulations; inasmuch as orbital decay is criti- 
cal for many of these systems, the use of rigid potentials may 
be problematic in some cases. 

As the third column of Table 2 shows, a wide range of cri- 
teria were used to define an acceptable match to the observa- 
tions. Many studies still seemed more focused on morphol- 
ogy than kinematics, and less than half presented compelling 
quantitative comparisons between models and data. Matches 
labeled "kin2" used 2-D kinematical information, and pre- 
sented the models and data in such a way that direct and un- 
ambiguous comparisons could easily be made - ideally, the 
model and data were overplotted, or at a minimum plotted 
to the same scale and orientation. Matches labeled "kinl" 
used 1-D data (e.g., long-slit spectra or mean velocities plot- 
ted as functions of a single coordinate), but again compared 
models and data directly. A few studies, labeled "gen" (for 
genetic), evaluated matches numerically; these will be dis- 
cussed in § 4.2.1. Most of the remaining studies, while draw- 
ing on spatially-resolved kinematic information, presented es- 
sentially qualitative comparisons between models and data; 
these matches are labeled "qual" in the table. This designa- 
tion is rather broad, ranging from studies which plotted mod- 
els and data on different scales to studies which matched gen- 
eral kinematic trends or asserted, without providing quanti- 
tative evidence, that the model matched the data. Finally, a 
few studies which matched morphology only are designated 
"morph"; these are included when they served as precursors 
to more comprehensive modeling efforts. 

The present study, while restricted to artificial data, closely 
parallels the approach of Hibbard & Mihos (1995). We depart 
from them in treating the pericentric separation p as a free pa- 
rameter, and in not using fully self-consistent simulations to 
finalize the models. The latter, of course, is deliberate; one 
of our goals was to see if an approach combining test-particle 
disks with self-consistent halos can recover the encounter pa- 
rameters of interacting disk galaxies. In practice, we envision 
using Identikit models to jump-start fully self-consistent sim- 
ulations. We strongly concur that kinematic data provide an 
acid test which any dynamical model must pass (Toomre & 
Toomre 1972; Borne & Richstone 1991), and that large-scale 
tidal features are the key to unlocking the dynamical history of 
galactic encounters. Finally, we emphasize that direct and un- 
ambiguous comparison between the simulations and the "ob- 
servational" data was a necessary ingredient of our approach. 
Overplotting the particles on the data, as in Fig. 1, is an effec- 
tive way to present such comparisons. 

4.2. 1 . Genetic Algorithms 

To date, most attempts to model interacting galaxies have 
relied on expert judgement in selecting initial conditions and 
identifying good matches between simulations and observa- 
tions. Recognizing the considerable labour involved, several 
groups have tried to automate the modeling process (Wahde 
1998; Theis & Kohle 2001; Gomez et al. 2002). The pro- 
posed algorithms have two essential components. First, they 
must replace the subjective comparison of the simulation par- 
ticles {3^) and the observed data (^) with an objective crite- 
rion ,^{,'^,&) measuring goodness of fit. Second, they must 
perform an efficient search of a very large parameter space. In 



view of the number of parameters involved, a blind search is 
impractical; these groups have adopted strategies mimicking 
biological evolution, generally known as genetic algorithms 
(Holland 1975). 

Genetic algorithms create a population of A^pop individuals, 
each representing a possible solution to the problem at hand; 
in this case, an individual defines a set of initial conditions 
and viewing parameters. The evolutionary fitness of individ- 
ual i is determined by using its initial conditions and viewing 
parameters to produce a particle distribution which is eval- 
uated using Once all A^pop individuals have been 
evaluated, the fittest among them are bred together to form 
a new generation, and the entire process is repeated. After 
A'gen generations, the population converges toward a nearly- 
optimal ensemble, with the fittest individual representing the 
best approximation to the desired solution. 

Genetic algorithms for modeling interacting galaxies have 
been tested on artificial data (Wahde 1998; Theis & Kohle 
2001; Gomez et al. 2002) and applied to real data for 
NGC 4449 (Theis & Kohle 2001) and NGC 5194/95 (Wahde 
& Donner 2001; Theis & Spinneker 2003). Typical values 
of A^pop ~ 10^ and A'gen ~ 10^ imply that ~ 10"* individuals 
must be evaluated to obtain a good match; with test-particle 
methods, this can be done in a few hours of CPU time. The 
Identikit methodology (§ 2) could be combined with a genetic 
algorithm, improving the treatment of orbital decay and sub- 
stantially reducing the CPU time required to find a match. 

However, the output of a genetic algorithm will be no bet- 
ter than the evaluation function ^(^, ^) used to determine 
fitness. The simplest approach is to coarsely grid ^ on the 
plane, and compare the result with an equally coarse 
gridding of 3i. More recent implementations incorporate ve- 
locity information as well, and there's no reason why grid- 
ding can't be extended to 3-D (X, F, V) data. But the first two 
limitations of the Identikit method noted above also apply to 
existing evaluation functions. Features due to self-gravity in 
real systems can be discounted by an expert when attempt- 
ing to fit a test-particle model, but may mislead an objective 
evaluation function, lowering the fitness of accurate solutions. 
And while HI is a good tracer of kinematics, its distribution in 
tidal features is often quite irregular; for example, a tail may 
appear as a series of clumps rather than a connected structure. 
An expert can recognize such tails as connected structures, but 
an automatic procedure may reject solutions which populate 
them with a smooth distribution of particles. To address these 
problems, recent genetic algorithm implementations include 
routines for masking or weighting the observational data; it 
remains to be seen if these techniques make genetic algo- 
rithms competitive with human experts. 

4.3. Are Models Unique? 

Can a dynamical model reproducing the morphology and 
kinematics of an interacting pair of galaxies be considered 
unique? Claims to this effect occasionally appear (e.g.. Borne 
1988a,b; Theis & Kohle 2001); skeptics, paraphrasing John 
von Neumann^, may be tempted to reply "with that many pa- 
rameters I could fit an elephant". The 16 parameters intro- 
duced in § 1 are all physically motivated and necessary to de- 
scribe a collision of two disk galaxies in 3-D; if a large num- 

See Dyson (2004) for one version: "with four parameters I can fit an 
elephant and with five I can make him wiggle his trunlc". John von Neumann 
could presumably fit a whole herd of elephants with the parameters used to 
describe a single galactic encounter! 
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ber of parameters per se was really a flaw, models of spectral 
line formation in stellar atmospheres, requiring up to ^ 90 
abundance parameters, would be on shaky ground indeed! Yet 
claims of uniqueness seem overconfident. A particular match 
may be unique within the universe of possibilities defined by a 
given model, implying that all the parameters appearing in the 
model can be determined within reasonable accuracy, and that 
no set of parameters outside this tolerance range yield as good 
a match. However, this is not the same thing as uniquely de- 
termining the dynamical state of a pair of colliding galaxies, 
which is specified by the distribution function /(r, v). 

The problem of dynamically modeling isolated early-type 
galaxies, which has generated an extensive literature, illus- 
trates some of the difficulties involved in determining /(r, v). 
Recent studies (e.g., van de Ven et al. 2008; van den Bosch 
et al. 2008) use Schwarzschild's (1979) metiiod to fit mod- 
els with equilibrium distribution functions depending on three 
integrals of motion to stellar velocity data obtained from in- 
tegral field spectroscopy. Such models are quite successful 
at describing the orbital structure of galaxies and diagnosing 
the presence of black holes and dark halos. But not many of 
these models are truly unique; the orientations and intrinsic 
shapes of axisymmetric models appear uncertain (Krajnovic 
et al. 2005; van den Bosch et al. 2008), while triaxial models 
pose additional difficulties (van de Ven et al. 2008). Deter- 
mining /(r, v) for one galaxy is hard; doing so for a pair of 
galaxies seems harder still. 

In the context of the model described here, it may not even 
be obvious that we can hope to constrain even the most basic 
parameters. The orientation of a single disk, unless it hap- 
pens to be exactly edge-on, cannot be determined from a data 
cube F(X,F,y); for example, a disk tilted by 45° to the line 
of sight produces exactly the same data cube as one tilted by 
135°. (This ambiguity can be resolved by using dust lanes to 
determine which side of the disk is closer, although doing so 
goes beyond a strictly kinematic approach to galaxy model- 
ing; see, e.g., the discussion in Struck et al. 2005.) By ex- 
tension, a pair of disks yield a four-fold degenerate solution, 
since each disk has two possible orientations. However, tidal 
interactions between disks break this degeneracy; as Fig. 14 
shows, two disks which initially display identical F(Z,F,y) 
distributions can be differentiated after an encounter. In this 
example, as in most of our Identikit solutions, the disk orien- 
tations (/i,a)i) and (i'ljCOi) are strongly constrained; a little 
"wiggle room" exists around the actual values, but no radi- 
cally different choice can reproduce the "observed" morphol- 
ogy and kinematics. 

Taken together the results in § 3.3 show that almost all of 
our 18 "good" and 12 "fair" fits accurately constrain the disk 
orientations, viewing direction, time since pericenter, pericen- 
tric separation, and scale factors. We did not examine errors in 
center-of-mass position and velocity, but the offsets in our fits 
were always very small, so presumably these parameters are 
also well-constrained. While some solutions are better than 
others, all 30 of these fits appear to be unique; the estimated 
parameters are always in the neighborhood of the actual val- 
ues, and often very close indeed. 

The 6 "poor" fits tell a different story. Object 16 (see 
Fig. 2), a wide retrograde encounter "observed" long after 
pericenter, has such attenuated tidal features that fitting this 
system is much like fitting a pair of disks before their first 
encounter; our fit is clearly degenerate. As § 3.5 describes, 
object 6 also produced non-unique fits. Other poor solutions 
may also be degenerate, in view of the large misaUgnments 



in viewing direction and disk orientation they display (see 
Figs. 6 and 7). 

The fact that we were able to produce well-constrained and 
reasonably accurate solutions for 30 out of 36 systems de- 
pends critically on the quality of the "observational" data we 
used. Data with lower resolution and added noise would yield 
a less favorable outcome; in particular, incomplete velocity 
information can seriously compromise a solution. Models of 
NGC 7252 illustrate this point; early attempts using very hm- 
ited velocity data concluded that this object resulted from a 
merger of two retrograde disks (Borne & Richstone 1991), 
but once high-quahty HI data was available a very convinc- 
ing solution involving a collision of two direct disks was soon 
found (Hibbard & Mihos 1995). Fortunately, instruments like 
the Expanded Very Large Array now coming on-line should 
make it far easier to obtain the detailed velocity information 
needed to accurately constrain dynamical models of interact- 
ing galaxies. 

However, even a surfeit of data will not guarantee an un- 
ambiguous model in every case. One example may be the 
"Whirlpool" galaxy, NGC 5194/95; after three and a half 
decades of modeling based on ever-better HI data there is still 
no consensus regarding the number of passages required to 
account for NGC 5194's optical morphology and extended 
HI tail (e.g. Durrell et al. 2003, and references therein). Mod- 
els in which NGC 5195 makes only one passage can account 
for many of the observations but have trouble matching the 
tail velocities, while models in which this galaxy makes two 
passages match the velocities better but yield double tails (see 
Salo & Laurikainen 2000, Fig. 4). One possible solution to 
this puzzle is to assume that NGC 5 194's HI disk was warped 
even before its first and only encounter with NGC 5195; this 
could tilt the tail with respect to the inner disk in such a way 
as to match the observed velocities. Of course, pre-existing 
warps vastly complicate questions of uniqueness! In the spe- 
cific case of NGC 5194/95, stronger evidence in favor of a 
single passage seems necessary to justify introducing such a 
warp; better data on stellar velocities may help (Durrell et al. 
2003). 

Our limited knowledge of dark matter is a more fundamen- 
tal barrier to the creation of truly unique models for inter- 
acting galaxies. Models of isolated galaxies in equilibrium 
can represent the dark matter as a potential imposed by hand. 
In contrast, models of interacting systems should really treat 
dark matter as a full participant in the dynamical equations^. 
Given how little we know about dark matter, a good deal of 
guess-work will be part of any such model, and this is un- 
likely to change unless we can somehow measure the detailed 
distribution and kinematics of dark matter in individual galax- 
ies. Nonetheless, models of interacting galaxies have already 
set limits on the radial structure of dark halos (Dubinski et 
al. 1996; Mihos et al. 1998; Barnes 1999; Springel & White 
1999), and further studies may yield additional insights into 
the nature, distribution, and dynamics of dark matter. 

While it may never be possible to determine a unique distri- 
bution function /(r, v) for a specific pair of interacting galax- 
ies, there are still compelUng reasons to construct detailed 
models. A good model of a galactic collision can serve as 
a unifying hypothesis, providing a context to integrate and in- 
terpret a variety of observations; for example, a time-line of 

* To be sure, not every study in Table 2 did this, but the speed of mod- 
em computers leaves little excuse for not treating the dark matter self- 
consistently in future work! 
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Fig. 14. — Identikit simulation showing how a tidal encounter breaks degeneracy in disk orientations. The system is viewed along the orbital axis. The disks 
have inclinations i'l = 45° and 12 = 135° and arguments a)i = 0)2 = 30°. Left: one time unit before pericenter; note that both disks present the same appearance 
on the {X,Y), (V,Y), and {X,V) planes (top left, top right, and bottom left, respectively), although the (X,Z) projection (bottom right) shows that they do not 
have the same orientation in space. Right: one time unit after pericenter; the difference between the two disks is now evident in all four projections. 



past encounters may help to understand interaction-induced 
star formation. Conversely, a bad model which contradicts 
other lines of evidence can be rejected. Viewed as a means 
to an end, dynamical modeling stands to teach us a good deal 
about galactic encounters and galaxies in general. 

4.4. Extensions 

While the current Identikit system is essentially an interac- 
tive modeling tool, the basic approach offers other possibil- 
ities. For example, once the orbital parameters p, e, and jj., 
time t, viewing angles {9x,9y,9z), and scale factors .if and 
'f have been fixed, it is possible to invert the mapping from 
initial disk orientation (/, (o) to the data cube F{X,Y,Z), and 
ask which (;, co) values allow a disk to populate a given point 
in the {X,Y,V) space. As a rule, each {X,Y,V) point cor- 
responds to some range of (o) values, but different points 
along one tidal structure should all be populated by the same 
disk, so this range can be constrained. This approach could 
effectively automate the process of finding disk orientations 
once other parameters have been selected; it could particu- 
larly useful in establishing confidence limits for key parame- 
ters once an initial match has been determined. 

The Identikit technique has some interesting similarities to 
Schwarzschild's (1979) method for building triaxial equilib- 
rium systems. In a nutshell, Schwarzschild started with a 
stationary mass model, populated it with all possible orbits, 



and figured out how to add up the time-averaged density dis- 
tributions of these orbits so as to recover the original mass 
model. We start with a time-dependent mass model generated 
by an encounter of two spherical systems; each sphere is ini- 
tially populated with all possible circular orbits. But instead 
of finding a weighted sum of these orbits which reproduces 
the mass distribution, we try to select one co-planar family of 
orbits from each sphere to match the observed kinematics and 
morphology of interacting galaxies. Perhaps some variant of 
the algorithms used in Schwarzschild's method to determine 
orbital weights could be applied to the problem of modeling 
interacting galaxies. 

We hope to implement some of these extensions in Iden- 
tikit 2. 
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APPENDIX 
A. IDENTIKIT MODEL CONSTRUCTION 

To set up initial conditions for Identikit models we need to construct spherical equilibrium N-body models with mass profiles 
m{r) given by (2). We use Eddington's (1916) formula (e.g., Biimey & Tremaine 1987, p. 236) to compute the distribution 
function f{E). Kazantzidis et al. (2004) used this approach to construct equilibrium halos with density profiles described by a 
three-parameter family of models (e.g., Zhao 1996). We go one step further by correcting for the finite resolution (i.e., "soften- 
ing") of the N-body force calculation. 
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We account for the effect of Plummer softening (Aarseth 1963; Hernquist & Barnes 1990) in the N-body simulations by 
introducing a quasi-empirical transformation of the total mass profile (Barnes, in preparation): 



m{r) 



l + (2/3)('^/«)(e/r)'^ 



(a/K) 



m(r) . (Al) 



where a is the logarithmic derivative of the density profile as r ^ 0, the parameter e is comparable to the softening length e, and 
the parameter k adjusts the shape of the transition near r e. This smoothed mass profile is then used to compute the potential: 

^=G^, (A2) 

where <I> ^ as r ^ oo. After expressing the original density profile p(r) = {AKr^)^^dm/dr as a function of this potential 4>, we 
compute the distribution function: 



f(E) = ^ / d^((i>-E)-^'^— (A3) 

Once f{E) has been calculated, generating an N-body reaUzation with A^sphr massive particles is straightforward. For each 
particle ;, let f; and v, be two vectors drawn from a uniform distribution on the unit sphere S^. Then the position of particle / 
is r, = r/f;, where r, is chosen by drawing a random number x from a uniform distribution in the range [0, w(oo)] and solving 
m{ri) = X, and the velocity is v,- = v,- v,-, where v,- is chosen from the speed distribution v^f{^v^ + '&(?",)) using rejection sampling 
(von Neumann 1951). The particle masses are w; = m{°o) /A^sphr- 

The Attest test particles are initially placed in a single disk. If the test particle distribution follows md(r), the initial radius qi 
of particle i may be selected by drawing a random number x from a uniform distribution in the range [0,OT(j(°°)] and solving 
mi{qi) = X. However, we prefer to bias the distribution by a factor of to improve disk sampUng at large r (see § 2). Let 

n{r)= I dxX^KiX), (A4) 
Jo 

where m'^{r) = dm^/dr, then qt is selected by drawing x from [0, TJ (oo)] and solving TJ (^,) = x. The orbital velocity v,- of particle 
i is calculated using the smoothed profile: 

Vi = ^J Grn{qi) / qi . (A5) 

Note that test particles are placed on exactly circular orbits, creating a perfectly "cold" disk. Finally, the position r, and velocity 
V, of particle i are rotated to aUgn r, x v, with a normalized angular momentum s, drawn from a uniform distribution on S^. 

B. RANDOM MERGER MODELS 

Our galaxy construction procedure has some elements in common with McMillian & Dehnen (2007). Like them, we compute 
isotropic distribution functions fh{E) and fk{E) for the bulge and halo, respectively, by approximating the disk's gravitational 
field with its spherically averaged equivalent. Unlike them, we use the resulting bulge and halo "as is", without first adiabatically 
imposing a flattened disk potential; the response of the bulge and halo to such an adiabatic transformation is so subtle that a 
good approximation to equiUbrium is possible without it. This makes our procedure quite fast. What follows is a brief technical 
description of our procedure; a full discussion and numerical tests will be presented elsewhere (Barnes, in preparation). 

The bulge follows a Hernquist (1990a) model out to a radius b\,, and tapers at larger radii to avoid placing a small number of 
particles at extremely large distances: 

«bWb 1 . , 

r, r<bb 



In ria-^ + r)- 

2 (Bl) 



Pb{r) = < 

where m\, is the bulge mass, and is fixed by requiring that Pb{r) be continuous at r = fob- For b\,^a\„ the slope dph/dr is also 
continuous at r = foj,- 

The halo follows a Navarro et al. (1996) model out to a radius foh. and tapers at larger radii as proposed by Springel & White 
(1999): 

»^h(ah) 1^ 

4;r(hi(2)-i)K-+ah)2' " ' 
Ph{r) = { p (B2) 



Ph 

where OTh(<3!h) is the halo mass within radius ah, and and j3 are fixed by requiring that Ph('') and dp^/dr are both continuous at 
r = foh- The halo is tapered more abruptly than the bulge to tame the logarithmic divergence of the standard Navarro et al. mass 
profile as r ^ oo. 

The disk has an exponential radial profile (de Vaucouleurs 1959a,b; Freeman 1970) and a sech^ vertical profile (van der Kruit 
& Searle 1981): 

Pi{q,(^,z) = -^.-^/"^sech^ (^) , (B3) 
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where {q = ^/x^ +y^,^,z) are cylindrical coordinates and is the total mass of the disk. 
For the bulge and halo, cumulative mass profiles are obtained by integrating (Bl) and (B2): 

rr rr 

mb{r)= dx4nx^Pb{x), and mi,{r) = dx^nx^ ph{x) ■ (B4) 
Jo Jo 

We use the cumulative mass profile for an infinitely thin disk, 

/Wd(»-)=/Wd(l-e-''/''d)(l+r/ad), (B5) 

which is adequate for our purposes since mb(r) » m^ir) at small r. These functions are summed to get the total mass profile 
m{r), the smoothed profile m(r) is computed using (Al), and the potential 4>(r) is computed using (A2). We then express p\, and 
Ph as functions of <I>, insert these functions in (A3) to obtain fb{E) and fh{E), and construct N-body reahzations of the bulge and 
halo following the procedure described in Appendix A. 
The disk is reaUzed by sampling an approximate distribution function 



fd{q,(l>,Z,Vq,V^,Vz) °= Pd{q,(j),z) Jf 



a<i,{q) 



(B6) 



where Vq, and Vz are velocities in the radial, azimuthal, and vertical directions, respectively. The function v{q) is the mean 
rotation velocity, while C7q(^), (q), and (7z{q) are dispersions in the radial, azimuthal, and vertical directions, respectively. The 
function ^{x) is a Gaussian, while Jf{x) resembles ^{x) but cuts off faster for large \x\: 



^(jc)oce-2^, ^(jc)oce-2(Vc)2-|(Vc)^, (B7) 

where c is fixed by requiring / dxx^ J^{x) = f dxJf{x). This function is used instead of a Gaussian to avoid overpopulating the 
high- velocity tail of the distribution. 
The local circular velocity Vc{q) is given by 

vl{q)=G^^+q-^, (B8) 
q dq 

where ms(<?) is the smoothed spheroid (bulge + halo) mass profile and <I>d is the potential due to the disk. To compute ms(r) we 
insert the spheroid profile ms{r) = m],{r) + mh{r) in (Al). Our expression for <t>d explicitly takes "softening" into account: 

d<i>i Grrid ke-''^<>Ji{kq) 

— dk -—^ T-— , (B9) 



dq al Jo (a, 2 + ^2)3/2 



where Ji (x) is the cylindrical Bessel function of order one, and setting ei = ^^Je^+z^ allows - in an approximate way - for the 
finite thickness of the disk. 

The vertical dispersion is given by the solution for an isothermal sheet (e.g., Binney & Tremaine 1987, p. 282): 

az{q) = ^KGzd^q), (BIO) 

where T,{q) = f dzpd{q,z) is the surface density of the disk at cyUndrical radius q. The radial dispersion is then determined by 
fixing the ratio (Jq/Gz- 

Cyq(?)=M(^)cTz(?)= + ) Oziq), (BID 

where is a scale parameter comparable to ad- The function ^i{q) is chosen to make Oq/Oz ~ 2 in the body of the disk 
- roughly matching the solar neighborhood value (e.g., Dehnen & Birmey 1998) - while letting Cq/ Cz — > 1 for g — > 0. The 
azimuthal dispersion is related to the radial dispersion (Birmey & Tremaine 1987, p. 203): 

M^) = ^^^^i^)^ (B12) 



where Q.{q) = Vc{q)/q is the circular orbital frequency and K{q) = y^HM-^^i?/^ is the epicychc frequency. 

Finally, the mean rotation velocity v{q) is determined using the axisymmetric Jeans equation (e.g., Binney & Tremaine 1987, 
p. 198): 

V'{q)=vl{q) + a'^{q) (^l -^^ - aj{q) + a^{q)q^ (B13) 

The parameter values needed to completely define the galaxy model are: 

OTb = 0.0625, ab = 0.02, fob =4.0, 

md = 0.1875, ad=l/12, Zd = 0.0075, 

mh(ah)=0.16, ah = 0.25, fob = 0.98015, (B14) 

£ = 0.0075, e = 0.01 15, es = 0.0115, 

^CT = 0.075, k:= 1.975, k:s = 2.025. 
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Fig. 15. — Circular velocity profiles for the disk galaxy model. Left: circular velocities for the bulge (red), disk (grey), and halo (blue). The total circular 
velocity is also shown (black). Right: total circular velocity (black), mean rotation velocity (red), and circular velocity for the spherical mass model (blue). 



A few remarks about these parameters are in order First, tapering the bulge as in (Bl) with b\, = 200flb reduces the total bulge 
mass by ^ 0.5%; to correct this, the value of mb actually used in (Bl) is adjusted upward accordingly. Second, the primary halo 
mass parameter is mh(flh); the halo taper radius bh is adjusted to make the total halo mass mh{°°) — = 1-0. Third, the softening 
parameter actually used in the N-body calculations is e; the values of e and K listed here are chosen by computing <I>(r) using 
(A2) and comparing the result to an N-body calculation. Fourth, the parameters Eg and fCj used to compute the smoothed spheroid 
mass profile Ms(r) are likewise chosen by comparison with an N-body calculation. 

Fig. 15 presents circular velocity profiles for the galaxy model adopted here. The left-hand panel shows profiles for the 
individual components, computed taking softening into account as described above. Also shown is the total circular velocity 
Vc{q) given by (B8). The right-hand panel again shows Vc(^) and compares it with the mean rotation velocity v(^) given by (B13) 
and the circular velocity for the equivalent spherical mass model given by (A5). 

To check this model, we constructed a realization with A^b = 16384 bulge particles, = 49152 disk particles, and A^h = 65536 
halo particles. This system was then evolved in isolation for 10 time units, using a hierarchical N-body code^ with an accuracy 
parameter = 1, quadrupole-moment corrections (Hernquist 1987), a Plummer softening length £ = 0.0075, and a leap-frog 
integrator with a time-step Af = 1 /256. During the first 0.25 time units the ratio of kinetic to potential energy, T /U, fell from an 
initial value of 0.4980 to 0.4945; it then fluctuated around this value with an amplitude of ^ 0.003. This initial drop indicates 
that the model was not started in perfect equilibrium, but the implied rearrangement of mass is only a little larger than the 1 / ^/N 
(~ 0.0028) fluctuations occurring in an N-body system with this A^. Apart from transient spiral structure, this model showed no 
significant features until the disk begins to develop a bar at time f ~ 4. It's hard to completely suppress a weak bar instability 
in galaxy models with relatively massive disks like the one used here; however, this instability has little effect on the merger 
simulations since the galaxies interact with each other long before they would develop bars in isolation. 

The 36 random merger simulations were run with the same number of particles per galaxy and N-body integration parameters 
used in the test just described. Energy was conserved to ^ 0.05% even in the most violent encounters. Particle positions and 
velocities were output every Afout =1/32 time units, providing a large data base which could be used to construct random samples 
like the one in Fig. 2. 

C. VELOCITY SCALE BIAS 

As noted at the end of § 3.3, Identikit estimates of the velocity scale Y are typically ^ 10% too high. The value of Y is usually 
determined toward the end of the matching process; after values have been selected for most other parameters, we adjust Y to 
obtain a good overall match between the particles and the grey-scale images in the {X,V) and {V,Y) planes. There are a number 
of factors which may influence the choice for y, including the velocity difference between the galaxies and the characteristic 
velocities of tidal features. However, a key feature is the velocity widths of the galaxies; we generally try to adjust Y so that the 
particles span the full range of velocities present in each galaxy. 

While matching the velocity widths of the individual galaxies seems reasonable, it appears to be the source of the bias in f^. 
Velocities in the test-particle Identikit disks are not identical to those in the self-gravitating disks used in the random mergers; the 
former are perfectly "cold", while the latter have non-zero velocity dispersions. In addition, as the right-hand panel in Fig. 15 
shows, the circular velocity profile of the Identikit model rises more rapidly, peaks at a smaller radius, and then falls slightly 
below either the circular (vc(^)) or mean {v{q)) profiles in the self-gravitating model. 

Fig. 16 shows an Identikit fit to a self-consistent disk galaxy model. This disk is tilted by 45° to the line of sight; in fitting 
this model, the orientation and length scale factor were first set to their actual values. We then adjusted Yfn by eye, stretching the 
particle distribution to match the apparent velocity width of the grey-scale images in the {X,V) and {V, Y) projections. The value 
of Yfit selected in this manner is 10.9% too high, much as in the Identikit matches to interacting galaxies. In hindsight, we note 
that the {X,V) projection shows a handful of points falling outside the rotation curve for large \X\; a fit giving priority to these 
points would have recovered a more accurate value for 

Compared to its Identikit analog, the larger velocity width of the self-consistent disk is a direct consequence of the random 
velocities imparted to individual disk particles. To test this, we constructed self-consistent disks with different dispersions and 

' See http://www.ifa.hawaii.edu/faculty/barnes/treecode/treeguide.html for a discussion of this code, which generalizes earlier modifications (Barnes 1990) of 
the original tree code (Barnes & Hut 1986). 
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Fig. 1 6. — Identikit matcii to a single disk galaxy. The velocity scale factor was adjusted by eye to match the width of the velocity distribution. 



velocity scales, and selected best-fit values of 'f&t as in Fig. 16. The ratio 'fkt/^tme has a one-to-one relationship with the disk's 
velocity dispersion; the larger the dispersion, the greater the overestimate of i^. It's likely that by constructing Identikit models 
with random velocities comparable to those present in the self-gravitating disks we could largely remove this bias. However, 
"cold" disks may actually be more appropriate when fitting kinematics observed in cold gas tracers (e.g., HI or CO), since these 
have smaller velocity dispersions than most stellar components. 
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TABLE 1 

Ratios of ht/true parameter values for all 36 

MODELS. 



Parameter 




Quart 


Median 


3"" Quart 


time since pericenter 


t 


0.86 


1.02 


1.14 


pericentric separation 


P 


0.79 


1.01 


1.25 


length scale factor 




0.86 


1.01 


1.17 


velocity scale factor 


r 


1.05 


1.10 


1.15 



physical time T 0.78 0.95 1.02 

physical separation P 0.87 0.97 1.13 
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TABLE 2 

Dynamical models of interacting disk galaxies^. 



Svstem'' 


Arp 


VV 


Tvne 


Match^ 


Reference 


VV784 




784 


ring 


moq^h 


Toomre (1978) 










kinl 


StructMarcell & Higdon (1993) 


VV347, Aip 119 


119 


347 


ring 


qual 


Hearn&Lamb (2001) 


NGC 520 


157 


231 


S+S 


kinl 


Stanford & Balcells (1991) 


3C 48 






s+s 


kinl 


Scharwachter et al. (2004) 


NGC 672 /IC 1727 




338 


S+S 


kin2 


Combes et al. (1980) 


NGC 1143/44 


118 


331 


E+S 


qual 


Lamb, Hearn, & Gao (1998) 


IC 1908, AM 0313-545 






S+d 


qual 


Mihos et al. (1993) 


NGC 2207 /IC 2163 






S+S 


qual 


Elmegreen et al. (1995) 










kin2 


Struck et al. (2005) 


W 785, AM 0644-741 




785 


ring 


kinl 


Antunes & Wallin (2007) 


NGC 2442 / AM 0738-692 






S+d 


qual 


Mihos & Bothun (1997) 


nHz4 






ring 


qual 


Lynds & Toomre (1976) 


NGC 2782 


215 




S+?'' 


qual 


Smith (1994) 


NGC 2992/93 


245 




S+S 


kin2 


Due et al. (2000) 


NGC 3031/77 






s+d 


qual 


Thomasson & Donner (1993) 


NGC 3031/34/77 






S+S+d 


qual 


BrouiUet et al. (1991) 










qual 


Yun (1999) 


AM 1003-435 






S+S 


qual 


Giinthardt at al. (2006) 


NGC 3395/96 


270 


246 


S+S 


kinl 


Clemens et al. (1999) 


NGC 3448 / UGC 6016 


205 




s+d 


qual 


Noreau & Kronberg (1986) 


NGC 4038/39 


244 


245 


S+S 


morph 


Toomre & Toomre (1972) 










kin2 


van derHulst (1979) 










qual 


Mahoney et al. (1987) 










qual 


Barnes (1988) 










qual 


Mihos etal. (1993) 










kin2 


Hibbard (2003) 


NGC 4254/92?"= 






S+S 


qual 


Due & Boumaud (2008) 


NGC 4435/38 


120 


188 


so+s 


moq^h 


Combes et al. (1988) 










kin2 


VoUmer et al. (2005) 


NGC 4449 / DDO 125 






S+d 


gen 


Theis & Kohle (2001) 


NGC 4654/39?^ 






s+s 


kin2 


Volhner (2003) 


NGC 4676 


242 


224 


s+s 


morph 


Toomre & Toomre (1972) 










qual 


Stockton (1974) 










kinl 


Mihos et al. (1993) 










qual 


Gilbert & Sellwood (1994) 










kinl 


Sotnikova & Reshetnikov (1998) 










kin2 


Barnes (2004) 


NGC 5194/95 


085 


001 


s+so 


morph 


Toomre & Toomre (1972) 










kin2 


Toomre (1978) 










qual 


Hemquist (1990b), 










kin2 


Salo & Laurikainen (2000) 










gen 


Wahde & Donner (2001) 










gen 


Theis & Spinneker (2003) 










kin2 


Durrell et al. (2003) 


NGC 5216/18 


104 


033 


E+S 


qual 


Cullen et al. (2007) 


NGC 5394/95 


084 


048 


S+S 


qual 


Kaufman etal. (1999) 


AM 2004-662 






E+d 


kinl 


Dfazetal. (2000) 


NGC 6872 /IC 4970 






S+SO 


qual 


Mihos etal. (1993) 










qual 


Horellou & Koribalski (2007) 


NGC 7252 


226 




S+S 


qual 


Borne & Richstone (1991) 










qual 


Mihos et al. (1993) 










kin2 


Hibbard & Mihos (1995) 










qual 


Mihos et al. (1998) 


NGC 7714/15 


284 


051 


ring 


qual 


Smith & Wallin (1992) 










kin2 


Struck & Smith (2003) 


NGC 7752/53 


086 


005 


S+d 


kin2 


Salo & Laurikamen (1993) 



^ This table attempts to survey and characterize dynamical models of interacting galaxies 
which make significant use of kinematic constraints. Different authors often use very 
different criteria when imposing kinematic constraints, and published descriptions are 
sometimes ambiguous. No warranty of completeness is expressed or implied. 

In this column, a slash separates components of a given system, while a comma sep- 
arates alternate names. NGC numbers are used when available; if all components have 
NGC numbers, the full nixmber is given for the first galaxy, and only the last two digits 
for the rest. 

Briefly, "kin2" matches are constrained by 2-D kinematic data, "kinl" matches are con- 
strained by 1-D kinematic data, '"gen" matches used genetic algorithms, "qual" matches 
reproduce qualitative kinematic features, and "morph" matches reproduce morphology 
(listed only as precursors). See text for details. 
^ Type of companion is ambiguous. 
^ Identity of companion is ambiguous. 



